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1. Introduction 

Inclusive quark production through a decay of a heavy virtual photon, Z boson or T is a process 
of importance for QCD as the theory of strong interactions. Perturbative QCD (pQCD) provides 
theoretically clean prediction for the process (see, e.g. [^, |]]). 

Combined with the precise determination of the Z-boson decay rate into hadrons at LEP []|] 
this has led to one of the most precise determinations of the strong coupling constant a s (Mz). An 
alternative and also very precise determination of a s (Mz) as derived from a s (M T ) has been re- 
cently obtained from the &((xf) prediction §\§ for the ratio R r = ^■(■z^i+t/Tv^) an( * tne experimental 
determinations of 7? T by ALEPH, CLEO and OPAL collaborations (see, e.g. [Q]). 

Note that while the &{oQ ) predictions for R t are complete this is not the case for R(s) and the 
Z-decay rate. The missing pieces are related to so-called singlet diagrams (see Fig. 1 below). Note 
that while the top quarks can not be produced in Z-decays due to kinematical reasons, the (axial) 
singlet diagrams containing internal quark loops are not power suppressed (unlike similar loops 
for the vector singlet (and non-singlet) diagrams. This remarkable phenomenon first shows up at 
order aj and was first established and fully investigated in works [||, [7|]. The full account of singlet 
diagrams at order q s 3 was performed in papers [||, |9|, [T^] (vector case) and in [ 11 , 12|, [TJ, [l0|] (axial 
case). 

In the present work we present the results of the calculations of the order a* axial singlet 
contributions for the decay rates of the Z-boson as well as the vector singlet contribution to the 
cross section for electron-positron annihilation into hadrons. Note that we will not dwell on any 
phenomenological applications of our calculations as they have been recently discussed in some 



detail in [14] 



2. Preliminaries 

The interaction of the Z boson to quarks is described (in the lowest order approximation 
in the weak coupling constant) by adding to the QCD Lagrangian an extra term of the form 

M z Z'Va* with J® = Y,iWiYa{gi — gfjsJWi being the neutral quark current. As a re- 

sult, the hadronic decay rate of the Z boson (r|) including all strong interaction corrections may be 
viewed as an incoherent sum of vector (T^) and axial (P|) contributions. By the optical theorem 
both contributions can be conveniently related to the correlators of vector and axial vector quark 
currents. The general definition for the latter reads: 



with j v i = Yi Yn Yi = V]i an d jZ ,■ = Wi YsYi = ^|u- The corresponding absorptive parts are 
defined as follows: 

Rl( A (s) = 12tt3 nl£(-s - is). (2.2) 
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(a) (b) (c) (d) 



Figure 1: Examples of non-singlet diagrams (a), (b), where the two Z vertices are connected by a fermion 
line, and of singlet diagrams (c),(d), where the diagram can be split by only cutting gluon lines. The imag- 
inary part of the non-singlet diagrams gives R V / A - NS 7 while the imaginary part of the singlet diagrams is 
denoted by R V I A < S . 

The Z decay rate T(Z -»■ hadrons) = F (R V (M^) + R A (Mj)), where T = G F M|/ (2<W2) and 
R V I A can be expressed in terms of Rjj A defined in eq. (|L2|), namely 

'J u 

Similarly, the inclusive cross-section reaction of the reaction e + e~ annihilation into hadrons 
through the photon is described by the current correlation function 

U^q) = J dxe^(0\T[ ;7(*)JT(0)]|0) = (-g, v q 2 + q,q v )U EM (-q 2 ), (2.4) 

with the hadronic EM current 

i;=Iw^/ and R(s) = nxm EM (-s-ie), 

f 

with q f being the EM charge of the quark /. As a result, we arrive to the following representation 
for the ratio R(s) valid in massless approximation (precise definitions of R NS and R vs will be given 
below) 

R(s)=£q i q j Rl j (s) = R NS (s) + fa,)* R vs (s). (2.5) 

ij i i 

As the Z-boson is much heavier than all known quarks but the top one, it is natural 1 to neglect 
all power suppressed light quark mass corrections when dealing with 0(ttf) contribution to T|. It 
is customary to split Rjj A into two contributions as described in Fig. [I] 

= 8fjR NS (M 2 ) + Ry A - s (Ml) , (2.6) 

with the delta function Sfj = 5jj if both flavours i and j are light and Sfj = if either i or/and j 
refer to the top quark. In the non-singlet diagrams there is no top quark present in the fermion loop 
connecting the two external currents, because these diagrams have no physical cut and therefore 
have no imaginary part contributing to R NS (s = M|). This, together with the assumed masslessness 

'Mass corrections to both vector and axial vector correlator due to other massive quarks are dominated by the 
bottom quark and can be classified by orders in m|/M| and a s . Up to ^{afmf/M^) and &{a^m l/M%) they can be 
found in as well terms of order a^m^/M^ (const + log m^/M|) and afm^/Mf (const + log m\/ Mj) t hat arise from 
axial vector singlet contributions. Terms of order a^m^/M^ and afrn^/M^ can be found in jl5| ] an d Jlp| | respectively. 
Corrections of order ajm^/mj and a^m^/mj from singlet and non-singlet terms are known from fl[ g, [Tf] and jlo| ] 
respectively. 
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of all quarks but top leads to the factorized form of the non-singlet term in eq. ( |2.6| ). Note that 
the internal top quark loops like in diagram (b) of Fig. 1 still contribute to R NS (s = Aff) if the 
strong coupling is defined for 6 flavours. However, it is well known that such contributions could 
be naturally described (up to power suppressed terms) by transition from the full nf = 6 QCD to 
the effective massless nf = 5 one (see, e.g. [^] and references therein). In fact, the same is true for 
the vector singlet term (0/j below is defined as 1 if either i or/and j refer to the top quark and in 
all other cases) 

R V if(M z ) = (1 - e^)R V ' S (M z ) + 0{M 2 z /Mf). 

The corresponding massless calculations of R NS in order have been recently finished [fi], 0|. 
In what follows we concentrate on the singlet terms R v,s and R A S . 



3. 75 -treatment 

As is well-known the treatment of 75 within dimensional regularization is a non-trivial prob- 
lem by itself (for an excellent review see [|T8|]). Following works [11, |12|] in all our calculations 
we employ, in fact, two different definitions of 75. First, for all non-singlet diagrams completely 
anticommuting naive 75 have been used. Second, for singlet diagrams we employ essentially the 't 
Hooft-Veltman definition ^\ 



(3.1) 



where the current YiYpYvYpYi * s assumed to be minimally renormalized. 

1 — 4a s /3 + 0{a 2 s ) on the rhs of (3J) is necessary [2C, 21] 



The finite normalization factor ^ 



for the current d3~l| ) to obey the usual (non-anomalous) Ward identities which in turn are crucial in 
renormalizing the Standard Model. 



In principle, one could (and even have to!) use one and the same definition (3. 1 ) also for non- 
singlet diagrams. This would result to much more complicated calculations due to significantly 
longer traces encountered. Fortunately, it is not necessary because the factor ^ is chosen in such 
a way to restore the anti-commutativity of the 75 (for a detailed discussion, see pi]]). 



4. Vector @(a A ) singlet term R 



v,s 



From purely technical point of view the calculation of the massless five-loop diagrams con- 



tributing to nj/ is not much different from those contributing to IT 7 '^. Using the same methods 
as described in [Jj], ^] we have obtained (below a s = a s (n)/n and ji. is the renormalization scale in 
the MS scheme) 
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R v > s (s) 
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5. Axial vector &{a A ) singlet term Rf'f 

Due to the obvious property if all 4 indexes refer to the massless quarks and the 

fact that g u A + g A = g c A + g A = 0, we can write the axial singlet part of the Z decay rate as follows: 3 



R A f = R%< s -2R? l ?+l$£. 



(5.1) 



All diagrams contributing to the first two terms of (5.1) contain at least one top quark loop. The 
third term receives contributions by both the completely massless diagrams and those with top 
quark loop (the latter start from order a 3 , an example is given by Fig. 1 (d)). 

As Mz € 2M,, one can use the effective theory methods to compute top-mass-dependent 
diagrams as a series in the ratio ^Jj-. The procedure was elaborated long ago and successfully 
employed (see works [22, 11, |l2|]) to get all ingredients of eq. ( |5~I1 ) at order a 3 at leading order 
in \/M t expansion (still keeping all power non-suppressed terms, including those which depends 
on ln(/i 2 /M 2 )). From purely technical point of view the evaluation at order a 4 involves absorptive 
parts of five-loop diagrams with massless propagators and, in addition, absorptive parts of four-loop 
diagrams combined with one-loop massive tadpoles, etc. down to one-loop massless diagrams 
together with four-loop massive tadpoles. The latter have been computed with the help of the 



Laporta's algorithm [23] implemented in Crusher [24]. The massive tadpoles with number of loops 
less or equal three have been independently recalculated with the help of the FORM program 
MATAD [f25|]. Our results for R? t ' S , Rff, and R AS read 
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(5.3) 



2 Obvious, thanks to the existence of the unitary SU (n/) symmetry in the flavour subspace of the first 1 . . .«/ = 5 
massless quarks. 



3t> 



Note that separate terms on the rhs of ( 5. 1 ) are not scale-invariant, while their sum is HI 1L 13]. 
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2 2 

Here = a s {pL)/% in the effective (topless) «/ = 5 QCD, = hi^-, ^u, = hij^r, and M, is 
the pole top quark mass. In addition, = £(n) is Riemann's zeta function and a n = Li„(l/2) = 

Er=iV(2n 



Finally, setting = Mz, we arrive at the following numerical form of (5.5) 
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6. Conclusion 

All our calculations have been performed on a SGI ALTIX 24-node IB -interconnected cluster 
of 8 -cores Xeon computers and on the HP XC4000 supercomputer of the federal state Baden- 



Wiirttemberg using parallel MPI-based [|26j] as well as thread-based [ |27| ] versions of FORM 



For evaluation of color factors we have used the FORM program COLOR [E9]. The diagrams have 



been generated with QGRAF []30|]. This work was supported by the Deutsche Forschungsgemein- 
schaft in the Sonderforschungsbereich/Transregio SFB/TR-9 "Computational Particle Physics" and 
by RFBR grants 11-02-01196 and 10-02-00525. 

We thank P. Marquard for his friendly help with the package Crusher. 
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